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Predicting the local flux of dark matter particles is vital for dark matter direct detection experi- 
ments. To date, such predictions have been based on simulations that model the dark matter alone. 
Here we include the influence of the baryonic matter for the first time. We use two different ap- 
proaches. Firstly, we use dark matter only simulations to estimate the expected merger history 
for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we 
use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, 
we find that a stellar/gas disc at high redshift (z ~ 1) causes merging satellites to be preferen- 
tially dragged towards the disc plane. This results in an accreted dark matter disc that contributes 
~ 0.25 — 1 times the non-rotating halo density at the solar position. An associated thick stellar 
disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars 
can be separated from those that formed in situ, future astronomical surveys will be able to infer 
the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an 
equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its 
low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, 
increases the likelihood of direct detection at low recoil energy, boosts the annual modulation sig- 
nal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP 
mass to be determined (see contribution from T. Bruch this volume). 



Identification of dark matter 2008 
August 18-22, 2008 
Stockholm, Sweden 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



A dark matter disc in the Milky Way 



J. I. Read 



1. Introduction 

The case for dark matter in the Universe is based on a wide range of observational data, from 
galaxy rotation curves and gravitational lensing, to the Cosmic Microwave Background Radiation 
[7]; [3]. Of the many plausible dark matter candidates in extensions to the Standard Model, Weakly 
Interacting Massive Particles (WIMPs) stand out as well-motivated and detectable [7], giving rise to 
many experiments designed to detect WIMPs in the lab. Predicting the flux of dark matter particles 
through the Earth is key to the success of such experiments, both to motivate detector design, and 
for the interpretation of any future signal [7]. 

Previous work has modelled the phase space density distribution of dark matter at solar neigh- 
bourhood using cosmological simulations that model the dark matter alone (see e.g. [13], [12]). 
Here we make the first attempt to include the baryonic matter - the stars and gas that make up the 
Milky Way. The Milky Way stellar disc presently dominates the mass interior to the solar radius 
and likely did so also in the early Universe at redshift z = 1, when the mean merger rate in a ACDM 
cosmology peaks [10]. The stellar and gas disc is important because it biases the accretion of satel- 
lites, causing them to be dragged towards the disc plane where they are torn apart by tides. The 
material from these accreted satellites settles into a thick disc of stars and dark matter [17]. In this 
work, we quantify the expected properties of this dark disc. Its implications for direct detection 
experiments and the capture of WIMPs in the Sun and Earth are presented in the contribution from 
T. Bruch, this volume and in [2]. 

We use two different approaches. In the first approach (§||), we use dark matter only simula- 
tions to estimate the expected merger history of a Milky Way mass galaxy, and then add a stellar 
disc to measure its effect. This work is presented in detail already in [10]. In the second approach 
C©, we use cosmological hydrodynamic simulations of the Milky Way to hunt for dark discs. Both 
approaches are complementary in quantifying the expected properties of the dark disc. The former 
allows us to specify precisely the properties of the Milky Way disc at high redshift; the latter is 
fully self-consistent. 

2. Approach #1: adding a stellar disc to cosmological dark matter simulations 

We used a cosmological dark matter only simulation already presented in [4] to estimate the fre- 
quency of satellite-disc encounters for a typical Milky Way galaxy; further details are given in 
[10]. From the simulation volume, we extracted four Milky Way sized halos at a mass resolution of 
m p = 5.7 x 1O 5 M . The subhalos inside each 'Milky Way' and at each redshift output were identi- 
fied using the algorithm in [5] and then traced back in time to their progenitor halos as detailed in 
[10]. 

From our sample of four Milky Way mass halos and assuming that mergers are isotropic 
(which we must do since we do not know how the disc should align with the dark matter halo), we 
find that a typical Milky Way sized halo will have 1(±1) subhalo merge within 6 < 20° of the disc 
plane with v max > 80km/s; 2(±1) with v max > 60; and 5(±1) with v max > 40km/s. Away from the 
disc plane there will be twice as may mergers at the same mass [10]. It is important to stress that 
these numbers come from the distribution of fully disrupted subhalos, not the surviving distribution 
that is significantly less damaging. 

We then estimated the effect of a stellar disc on these mergers by running isolated disc-merger 
simulations. We set up our Milky Way (MW) model (disc+halo system) by adiabatically growing a 
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Figure 1: (a): The accreted stars (red) and dark matter (blue) at the end of a simulation where the LMC 
satellite merged at 9 = 10° to the Milky Way stellar disc; the black contours show the underlying Milky 
Way stellar distribution, (b): The corresponding velocity distribution in (rotation velocity) at the solar 
neighbourhood; the underlying Milky Way dark matter halo is shown in green, (c): The dark matter disc 
to dark matter halo density ratio Pddisc/Phalo as a function of height above the disc plane (also for 8 < 
R < 9 kpc), for LMC merging at 9 = 10°, 20°, 40° and LLMC at 9 = 10° to the Milky Way disc plane, (d): 
The satellite-disc inclination angle 9 as a function of time; the lines are truncated when the satellite is fully 
accreted. 

disc inside a spherical halo, as detailed in [10]. We chose three models for our satellite, but present 
just two here: LMC with v max = 60km/s, and LLMC with v max = 80km/s; these were set up as 
scaled versions of our MW model. We chose a wide range of initial inclination angles to the disc 
from 6 = 10 — 60°, one retrograde orbit, and range of pericentres and apocentres. The simulations 
were evolved using the collisionless tree-code, PkdGRAV [14]. The final evolved systems were 
mass and momentum centred using the 'shrinking sphere' method described in [11], and rotated 
into their moment of inertia eigenframe with the z axis perpendicular to the disc. 

The results are shown in Figure [[]. The left panel shows the accreted stars (red) and dark 
matter (blue) at the end of a simulation where the LMC satellite merged at 6 = 10° to the disc. 
Both the stars and the dark matter have settled into accreted discs. The middle panel shows the 
corresponding velocity distribution in (rotation velocity) at the solar neighbourhood (a cylin- 
der 8 < R < 9 kpc, |z| < 0.35 kpc). The underlying dark matter halo is shown in green and is not 
rotating; the accreted stars and dark matter (red and blue) have kinematics similar to that of the un- 
derlying stellar disc (black). The right panel shows the dark matter disc to dark matter halo density 
ratio Pddisc/Phalo as a function of height above the disc plane for selected merger simulations, 
as marked. As the satellite impact angle 6 is increased, the satellite contributes less material to a 
dark disc. For 9 = 40°, the density at the solar neighbourhood is nearly flat with z and less than 
a tenth of the underlying halo density; there is correspondingly less rotation in this simulation. 
Summing over the expected number and mass of mergers, we find that the dark disc contributes 
~ 0.25 — 1 times the non-rotating halo density at the solar position [10]. It is important to stress 
that all satellites regardless of their initial inclination have some accreted material that is focused 
into the disc plane (see Figure [j]), right panel. As such, we expect that the accreted dark and stellar 
discs will comprise several accreted satellites; the most massive low-inclination mergers being the 
most important contributors. 

The accreted stellar disc shares similar kinematics to the dark disc. Depending on assumptions 
about the mass to light ratio of accreted satellites, these accreted stars can make up ~ 10 — 50% 
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( a) MW1 ( b) H204 (c) H258 ( d) H258dark 




Figure 2: (a-cj The distribution of rotational velocities at the solar neighbourhood for three simulated 
Milky Way mass galaxies MW1, H204 and H258. The lines show the dark matter (black), stars (red), a 
double Gaussian fit to the dark matter (blue; blue dotted), the material accreted from the four most massive 
disrupted satellites (green, blue, magenta, cyan), and the sum of all material accreted from these satellites 
(black dotted). The best fit double Gaussian parameters are marked in the top left, along with the redshift, 
Z. (d) As (a-c), but for the galaxy H258 simulated with dark matter alone, (e-h) The decay in radius r as 
a function of redshift z of the four most massive disrupting satellites in MW1, H204 H258 and H258dark. 
Where less than four lines are shown, these satellites accreted at redshift z > 3. The dotted sections show the 
evolution of the most bound satellite particle after the satellite has disrupted. 

of the Milky Way stellar thick disc [10]. (The lower end of this range is more likely given the 
observed properties of satellite galaxies in the Universe today.) If future surveys of our Galaxy can 
disentangle accreted stars in the Milky Way thick disc from those that formed in-situ, then we will 
be able to infer the velocity distribution function of the dark disc from these stars. 

3. Approach #2: cosmological hydrodynamic simulations 

We use three cosmological hydrodynamic simulations of Milky Way mass galaxies, two of which 
(MW1,H258) have already been presented in [9] and [6]. All three were run with the GASOLINE 
code [16] using the "blastwave feedback" described in [15]. MW1 had cosmological parameters 
(O m , £l A , a 8 ,h) = (0.3, 0.7, 0.9, 0.7) ; H204 and H258 used (Q. m , £1 A , a 8 , h) = (0.24, 0.76, 0.77, 0.73). 
At redshift z = the typical particle masses for dark matter stars and gas were: (Mj m ,M*,Mg as ) = 
(7.6,0.2,0.3) x 1O 5 M , with associated force softening: (£</,„,£*, £ gtu ) = (0.3, 0.3, 0.3) kpc. The 
analysis was performed as in §|2[ 

Figure || shows the distribution of rotational velocities at the solar neighbourhood (top panels), 
and the orbital decay of the four most massive satellites (bottom panels) for MW1, H204 and H258. 
The right most panels (d) and (h) show the results for the galaxy H258 simulated without any gas 
or stars - H258dark. 

Notice that, with the exception of H258dark, in all cases the dark matter requires a double 
Gaussian fit to its local velocity distribution. Once the baryons (the stars and gas) are included 
in the simulation, there is a local dark matter disc that lags the rotation of the thin stellar disc 
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by ~ 50 — 150km/s. The mass and rotation speed of the dark disc increase for the simulations 
that have more late mergers. MW1 has no significant mergers after redshift z = 2 and has a less 
significant dark disc, with rotation lag ~ 150km/s and dark disc to non-rotating dark halo density 
ratio Pddisc/Phalo = 0.23 (obtained from the double Gaussian fit). H204 and H258 both have 
extreme dark discs with Pddisc/Phalo > 1 and rotation lag ^ 60km/s; they both have massive 
mergers at redshift z < 1. 

4. Conclusions 

Low inclination massive satellite mergers are expected in a ACDM cosmology. These lead to the 
formation of thick accreted dark and stellar discs. We used two different approaches to estimate the 
importance of the dark disc. Firstly, we used dark matter only simulations to estimate the expected 
merger history for a Milky Way mass galaxy, and then added a thin stellar disc to measure its effect. 
Secondly, we used three cosmological hydrodynamic simulations of Milky Way mass galaxies. In 
both cases, we found that a typical Milky Way mass galaxy will have a dark disc that contributes 
~ 0.25 — 1 times the non-rotating halo density at the solar position. The dark disc has important 
implications for the direct detection of dark matter [2] and for the capture of WIMPs in the Sun 
and Earth (see contribution from T. Bruch, this volume). 
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